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Abstract 

It is generally expected that quantum gravity theory will bring the picture of a 
space-time foam at short distances leading to Lorentz Invariance Violation (LIV) 
manifested e.g. by energy dependent modification of standard relativistic dispersion 
relation. One direction of research, pursued intensively is to measure energy de- 
pendent time of arrival delays in photons emitted by astrophysical sources located 
at cosmological distances. This is tempered however by our ignorance of intrinsic 
emission delays in different energy channels. 

In this paper we discuss a test based on gravitational lensing. Monitoring time 
delays between images performed in different energy channels (e.g. optical - low en- 
ergy and TeV photons) may reveal extra delays due to distorted dispersion relation 
typical in LIV theories - a test which is free from the systematics inherent in other 
settings. 
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1 Introduction 



Despite tlie fact that quantum gravity theory still remains elusive, it is gen- 
erally expected that it will bring the picture of a space-time foam at short 
distances leading to Lorentz Invariance Violation (LIV) manifested e.g. by en- 
ergy dependent modification of standard relativistic dispersion relation |1|2] . 
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Several years ago it has been proposed to use astrophysical objects to look 
for energy dependent time of arrival delays [T]. Specifically gamma ray bursts 
(GRBs) being highly energetic events visible from cosmo logical distances are 
the most promising sources of constraining LIV theories [3 .4.5, 6j. Among other 
sources the BL Lac objects like Mk 501 are considered. It is this particular 
object from which 20 TeV photons were reported [2]. Such objects (also called 
blazars) have similar nature with quasars. In the papers quoted above some 
limits on LIV energy scale have been derived (and recently corrected [^). 

The idea of searching for time of flight delays is tempered however by our 
ignorance of intrinsic delay (at source frame) in different energy channels (see 
e.g. [1]). Recently it was also shown [S], that lack of detailed knowledge about 
cosmological model (in the context of accelerating expansion of the Universe) 
could be another source of systematic effects at high redshifts. 

In this paper we discuss a test based on gravitational lensing. Its idea has 
been noticed in [1] in the context of gamma ray bursts playing the role of 
high energy photons. Namely such a source located at a cosmological distance 
may undergo gravitational lensing by a galaxy lying closer to the observer 
along the line of sight with an encounter parameter small enough to produce 
multiple images p] (the so called strong lensing). Indeed all of known strong 
lensing systems [TOj have quasar as a source and a galaxy (in most cases 
elliptical) as a lens. The light signals emitted by a source will be seen by 
the observer delayed (achromatically in classical General Relativity) at the 
location of images. This opens up a possibility to study time delays induced 
by LIV. Namely, monitoring the time delays between lensed images performed 
in different energy channels (e.g. optical or gamma-ray — low energy and TeV 
— high energy photons) may reveal extra delays due to distorted dispersion 
relation typical in LIV theories. This test is free from the systematics inherent 
in other settings. Next sections will substantiate the argument further. 



2 LIV induced time delays in different cosmological models 

Following [2] let us consider a phenomenological approach for LIV by assuming 
the modified dispersion relation for photons in the form: 



where:e = ±1 is "sign parameter" [2j, ^„ is a dimensionless parameter. As 
a first guess one may assume Eqg equal to the Planck energy, = 1 and 
^2 = 10~^ [6]. The dispersion relation ([T]) essentially corresponds to the power- 
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law expansion (see [3]) so for practical purposes (due to smallness of expan- 
sion parameter E/Eqg) only the lowest terms of the expansion are relevant. 
Because in some LIV theories the odd power terms might be forbidden [TT] 
usually the cases of n = 1 and n = 2 are retained. It should be noted that in 
some string theories third order corrections appear as leading ones. 

The relation ([1]) leads to a hamiltonian 

n 
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from which time dependent group velocity v{t) = ^ can be inferred. 
The comoving distance travelled by photon to the Earth is 

^detection -2 l / 







where in the last equation a standard time-redshift parametrization was taken 
into account. H{z) here denotes the Universe's expansion rate (so called Hub- 
ble function). Starting from this point our considerations will have cosmologi- 
cal connotation. The reason for this is simple — the modifications due to LIV 
theories are really tiny, so one has to look for sources located at cosmologi- 
cal distances (such like quasars or gamma ray bursts) which are far enough 
to compensate for the smallness of LIV corrections. This means that cosmo- 
logical background geometry should be taken into account. From now on we 
will assume it to be a flat Friedmann-Robertson- Walker model with Vt^ = 0.3 
and f^A = 0.7 — the so called concordance model, as supported by observa- 
tions [T2]. Some alternatives to the concordance model could in principle be 
discussed as well [8]. 

Expressing group velocity in terms of redshift, we get 

t;(z) ~c(l + z)[l + e^^^^ f-^) (1 + ^)1 (4) 



Time of flight for the photon of energy E is equal to 



In the first term one easily recognizes the time of flight for photons in standard 
relativistic cosmology (i.e. without LIV). Due to very small magnitude of LIV 
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corrections it also fairly represents the time of flight for low energy photons. 
Therefore below we neglect LIV corrections at low energy. 

Consequently, the time delay between a low energy and a high energy photon 
is equal to 



where we restricted our attention to "infraluminal" motion of high energy 
photons (i.e. low energy photons arrive earlier to the observer). Generalization 
to "superluminal" motion is straightforward — the same value with just an 
opposite sign (time delays become early arrivals). 

The idea of observational strategy emerging form ([6]) is again simple: monitor 
appropriate (i.e. emitting both low and high energy photons) cosmological 
source at different energy channels and try to detect this time delay. Some 
attempts along this line have already been undertaken [Sfip^ . However there 
remains an indispensable uncertainty about intrinsic time delays: there is no 
reason for which low and high energy signal should be emitted simultaneously, 
and while detecting distinct signals (peaks in the light curve) at different 
energies we have no idea which one was sent first. 

Our method outlined below, invoking gravitational lensing allows to get rid of 
this ambiguity. Before describing it in subsequent paragraph let us recall now 
that in cosmology one distinguishes three types of distances: 

(i) comoving distance: 



by r{z) we denoted a reduced (dimensionless) comoving distance, i.e. a co- 
moving distance expressed as a fraction of the Hubble horizon = c/Hq, 

(ii) angular diameter distance: 




(6) 
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(iii) luminosity distance: 




(9) 
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Angular diameter distance is the one used in gravitational lensing theory (be- 
cause gravitational lensing deals with light deflection i.e. essentially with an- 
gles). The luminosity distance is a measure invoked while using standard can- 
dles (e.g. SNIa). The point is that both distance measures are related to the 
comoving distance by a 1 + 2; factor (see above). The comoving distance, then, 
is closely related to the time of flight r{z) = ct. It is in fact the distance to 
the source measured in light years. Therefore we can rewrite the time of flight 
in LIV theory ([5]) in terms of comoving distance: r^jvi^) = ctuv- Of course 
the comoving distance to the source is fixed — there are photons of different 
energies that travel with different speed. It is however useful (for later calcu- 
lations) to think as if they travelled with the same speed c but along different 
comoving distances r{z) and riiviz)- 



3 Gravitational lensing time delays 

Gravitational lensing of quasars and extragalactic radio sources at high red- 
shifts by foreground galaxies is now well established and has developed into a 
mature branch of both theoretical and observational astrophysics . Misalign- 
ment of the source, the lens and observer results typically in multiple images 
whose angular positions and magnification ratios allow reconstructing lensing 
mass distribution. In particular they provide an independent confirmation of 
dark matter in galaxies and became an important tool for investigating dark 
matter distribution. Another important ingredient of gravitational lensing is 
the time delay between lensed images of the source. This effect originates as a 
competition between Shapiro time delay from the gravitational field and the 
geometric delay due to bending the light rays and is best understood in terms 
of Fermat principle. In other words, the intervening mass between the source 
and the observer introduces an effective index of refraction, thereby increasing 
the light travel time. 

In general, the light travel time can be calculated as 
tCx") = ^ ^ \^(^ -Iff - ^( x") 

where: if and (3 are positions (as projected on the celestial sphere) of the 
image and the source, ipi^lf) is the projected gravitational potential (i.e. the 
actual potential integrated along line of sight), Di, Dg are angular diameter 
distances to the lens and the source located at redshifts zi and Zg respectively 
{Dis is the angular diameter distance between lens and source). From now 
on we adopt the notation (standard in gravitational lensing theory) where D 
denotes angular diameter distance, and subscripts refer to the components of 
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lensing system (i.e. the source, the lens or the observer). 



The lensing is called strong if source position happens to lie within the so 
called Einstein ring — the circle of a radius (defining the proper deflection 
scale of a given lens). In this case multiple images appear and since lensing 
galaxies are often ellipticals, the number of images is usually equal to four 
[TOfT^ (or five [15] - the issue of image multiplicity is discussed e.g. in [9]). 
However, the surprisingly realistic model of the lens potential is that of a sin- 
gular isothermal sphere (SIS) |9j. Indeed lensing by elipticals can be modelled 
by its variant called singular isothermal ellipsoid (SIE). Therefore for the pur- 
pose of illustrating our ideas we shall restrict our attention to the SIS model 
since generalization to SIE is rather straightforward and would not change our 
conclusions. 



The Einstein ring radius for the SIS model is: 

..^4.^^ (11, 



where a denotes one-dimensional velocity dispersion of stars in lensing galaxy. 
If the lensing is strong i.e. (3 := \ (3\ < then two co-linear images A and 
B form on the opposite side of the lens, at radial distances Ra = P + '&e and 
Rb = — P having time delays between the images: 



which according to the above mentioned relations for SIS model can also be 
written as 



2(1 + zi) DiD, ^ Stt^ a 
Ats^s = -^^^eP = j^nP- (13) 
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In the last equation f/ denotes the reduced comoving distance to the lens. The 
equation f|T2l) is commonly used by gravitational lensing community because it 
reduces time delay problem to relative astrometry of images, whereas (3 is much 
harder to asses (it is small in order for strong lensing to occur) and Einstein 
ring radius is not a directly observable quantity. However, the equation (fT3l) is 
more useful from the theoretical point of view. In particular it shows explicitly 
that the time delay between images is created at the lens location (f; factor). 
Let us stress again that this time delay is achromatic in General Relativity. 
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4 LIV induced time delays and gravitational lensing time delays 



Let us now imagine a source at cosmological distance emitting low energy and 
high energy (in TeV range) photons which undergoes gravitational lensing by 
a foreground galaxy. Let us also assume that LIV type distorted dispersion 
relation ([1]) holds in nature. The observer would notice again time delays 
between images, but this time it would be a combined effect of gravitational 
lensing and LIV. Therefore it would no longer be achromatic. This idea was 
formulated originally in [Ij but up to our best knowledge it has not been 
further developed. 

It is rather straightforward to calculate this by using the above mentioned 
fictitious "LIV comoving distance" rLiv{z), namely: 

^tuv,sis = -Tr'^Liv{zi)(3— 14 



where: 



n + 1 / E Y J n + zTdz' 



rLiv{zi)=ri + Ho^—[——] / ^, (15) 
2 XinEQc) J H[z') 



Because the LIV effect is extremely small, let us restrict further to the n = 1 



case: 



rw(.0=n + //o^y^^ (16) 



Now we can assume that observations in low energy would essentially provide 
time delay between images equal to Atsis, whereas monitoring of the same 
images in high energy (TeV) channel would provide AtLiv,sis- These two 
measurements would differ by 

AtLiV,SIS - Atsis = -jrf^^T^ / ;^ (17) 



H, c^EqgI H{z 







Rigorously one should rather calculate Atuy^sis{.Ei)—/S.tuv,sis{.E2) and write 
AE = El — E2 in the numerator of f[T7|) but due to our assumption that E2 
is many orders of magnitude smaller than Ei and thus having negligible LIV 
correction the above expression is a good appproximation. Let us make an 
estimate for the above LIV effect taking a real strong lensing system. HST 
14176+5226 can serve as an example. This system was discovered with the 
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Hubble Space Telescope [H] and further confirmed to be a gravitational lens 
[16]. The lensed source is a quasar at redshift Zg = 3.4 whereas the lens is an 
elliptical galaxy having redshift zi = 0.809. The lens model best fitted to the 
observed images was based on a singular isothermal ellipsoid [T7] giving the 
Einstein radius Be = 1".489 and /5 = 0".13 = 8.4 x 10"^ rad. 

Optical spectroscopy of the lensing galaxy in HST 14176+5226 system [18] 
provided measurements of the velocity dispersion in lensing galaxy. These mea- 
surements have been confirmed by Treu and Koopmans [19] who performed 
spectroscopic observations on Keck telescope as part of Lenses Structure and 
Dynamics (LSD) Survey. The result is a = 290 ± 8 km/ s. 

Substituting these data to (ITTl) gives ^tuv,sis — ^isis equal to 3.7 x 10^^ s 
for 5 TeV photons and 1.5 x 10~* s for 20 TeV ones. 

The model presented above was the simplest one just because its aim was 
to illustrate ideas. So was also the purpose of choosing the HST 14176+5226 
system as an example. In reality one would encounter systems with different 
numbers of images or different separations, e.g. SDSS J1004+4112 [20] with 
an image separation of 14". 6 or recently discovered system SDSS J1029+2623, 
where a quasar at Zg = 2.197 is doubly imaged by a massive galaxy cluster at 
zi = 0.55 with images separation of 22". 5 \21\ . 

However none of these effects is crucial to our general arguments because of 
their differential setting {At high energy vs. Atio^ energy )• High and low energy 
photons from each image travel along the same paths (respectively) thus suffer- 
ing the same shear effects (respectively). We do not mean here that one could 
disregard these effects but rather to strengthen the standard lore of lensing 
community to treat each lensing system separately in a detailed manner. 

Closing this section it would be interesting to ask how the LIV effects might 
modify image configurations. It could be suspected that they might do so since 
from the Fermat's principle perspective images are located at stationary points 
of the wavefront travel time functional (given by equation ffTOj) ). Therefore 
since LIV modifies time of flight in an energy dependent way (due to modi- 
fied dispersion relation) then one expects the images seen at different energies 
located at different positions. It is easy to see that for the SIS lens (general- 
izations to other mass profiles are also rather straightforward) the difference 
between Einstein radii for high and low energy photons A9e,liv '■= Oe,liv — Oe 
would be given by formula: 



where: l'^^\zi,Z2) := ^^mit^ • -^^^ reahstic lens configurations like HST 



A6e,liv — Ge 




) 



(18) 
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14176+5226 this would give negligibly small corrections of order 10~^^ arc sec. 
Hence even if LIV were operating this would not be able to change macro- 
images position in a detectable way. However it cannot be excluded that such 
minute differences could become relevant while studying caustic crossing possi- 
bly leading to different magnification patterns due to microlensing at different 
energies. 



5 Discussion and conclusions 

In this paper we discussed a method (first noticed in [T]) to test LIV effects 
by monitoring time delays between images of gravitationally lensed quasars 
in low and high energy channels. In standard theory (General Relativity) the 
result should be the same — gravitational lensing is essentially achromatic. 
On the other hand in the presence of LIV effects time delays loose this prop- 
erty — high energy photons should come at different times comparing with 
low energy ones. Therefore time delays between images should be different 
at different energies (e.g. optical or gamma-rays and TeV photons). We are 
tacitly assuming (following approach taken by the rest of LIV studying com- 
munity) that LIV effects are manifested only in high energy domain (where 
the small scale "foamy" structure of the space-time reveals itself) whereas 
the overall background geometry of space-time shaped by low energy con- 
tent of the Universe is that of General Relativity (more precisely — the flat 
Friedman-Robertson- Walker model as suggested by cosmological data) with 
light deflection (i.e. geodesic motions) defined in a standard way. 

Because this method is differential in nature, it gets rid of the assumptions 
about intrinsic time delays of signals at different energies. In fact time delays 
between images at different energies could be established in different exper- 
iments (at unrelated observing sessions) performed on given lensing system. 
The only demand is that they are accurate enough (done with a sufficient 
temporal resolution). Since the time delay between images is produced at the 
lens location, the result does not depend very strongly on the cosmological 
model. Lenses are located at modest redshifts where all realistic cosmological 
models essentially agree. 

One may ask if appropriate lensing systems (i.e. having sources emitting both 
low and high energy photons) exist. It is an observational fact that very high 
energy emission [E > 100 GeV) has been detected from over a dozen of blazars 
[22] which have similar nature with quasars. Quasars, on the other hand are the 
sources in all known strong lensing systems — CASTLES database contains a 
100 of such systems [10]. It is a matter of coordinating strong lensing surveys 
with experiments in high energy astrophysics (such like AGILE, GLAST or 
MAGIC experiments [23]) and the future will certainly bring the discovery of 
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lensed high energy source. Angular resolution of high energy experiments is 
gradually being improved. For example a recently launched AGILE instrument 
[21] has been designed to obtain accurate localization (~ 2' — 3' ) of transient 
events by the Gamma Ray Imaging Detector (GRID) - SA combination. So 
it is already close to typical image separation and definitely being able to 
see images like those from SDSS J1004+4112 or SDSS J1029+2623 as sepa- 
rated. Depending on exposure and the diffuse background its flux sensitivity 
threshold can reach values of (10 — 20) x 10"*^ photons cm~'^ at energies 
higher than 100 MeV with with an effective area above 200 om? at 30 MeV . 
Moreover it has excellent timing capability, with overall photon absolute time 
tagging of uncertainty below 2 /is and very small deadtimes (< 200 /is for the 
GRID, ~ 5 /is for the sum of the SA readout units, and ~ 20 /is for each of 
the individual Csl bars). In fact, AGILE instrument is optimized in the range 
below 1 GeV hence it is not representative to TeV range experiments needed 
to probe LIV but it clearly shows the the gradual improvement of sensitivity, 
timing and angular resolution in high energy astrophysics. 

An order of magnitude estimate for the effect discussed in this paper is not 
encouraging now. For a typical lensing system like HST 14176+5226 it is of 
order of nano-seconds. However, having in mind that high energy astrophysi- 
cal sources display rapid variability (indeed intrinsic variability in relativistic 
shocks powering these sources is enhanced by a Lorentz factor typically of 
order of 10^) and that e.g. light curves of gamma-ray bursts are already sam- 
pled with mili-second resolution (and AGILE went down to micro-seconds) 
one should not reject the idea presented above on the grounds that it is not 
within the scope of present-day observational technology. 
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